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I. INTRODUCTION 



Not only is the physical interpretation of the radiation of the acoustic waves in 
more than one dimension complex, but the effects of varying frequency and shading 
are hard to visualize when considering directionality and frequency properties. As 
Wood implied in his book, A Textbook of Sound [Ref. 1: p.449], the directional factor 
is very important in many acoustic applications such as speech, hearing and array 
design. 

The simplest array one can consider is the beam pattern of a dipole radiating in two 
dimensions. References to the case of the simple dipole can be found in virtually every 
textbook in acoustics. For a fairly recent example, see Malecki's book [Ref. 2: p.199]. 
Currently, the mathematical modeling of the simple acoustic arrays has been enhanced 
with the advent of the sophisticated computer software. Specifically, Ziomek 
[Ref. 3: p. 1 38] and, Embleton and Thiessen [Ref. 4: p. 1 124] have investigated N point 
source arrays, considering both amplitude and phase weighting. Although the 
mathematics and computer hardware exist, the visualization of the three-dimensional 
aspects of these arrays is difficult. By use of the computer graphics software the ability 
to interpret the significance of shading and phasing can be fully developed [Refs. 5,6]. 

To date, visualization of the three-dimensional beam pattern has been needed to 
understand the directional factor. Numerical computer output was essential to draw 
three-dimensional graphs. By changing two parameters (in three-dimensions, those 
could be x and y or 0 and (p, etc.), different generated values of the beam pattern, F, 
are obtained. 

One computer program to visualize the beam pattern in three-dimensions (P3D) and 
four computer programs to read values in detail (2DICIR, 2DPCIR, 2DIFIX, 
2DPFIX) were written. These computer models predict the directional factor and 
overall beam pattern with or without differences in the source strengths, inter-source 
distances, and phases of the elements of the arrays. 

Typical geometry for the beam pattern prediction is shown in Figure 2.3 on page 

16. 



The following definitions apply and will be used throughout: 
1 = source strength of the first source (source 1) 
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A = source strength of source 2 

B = source strength of source 3 

C = source strength of source 4 

D 1 = distance from the geometrical center to source 1 

D-, = distance from the geometrical center to source 2 

D 3 = distance from the geometrical center to source 3 

D 4 = distance from the geometrical center to source 4 

<p 2 = phase difference between source 2 and source 1 

(p 3 = phase difference between source 3 and source 1 

(p 4 = phase difference between source 4 and source 1 

8 = angle measured from the z axis toward the x-y plane 

<P = angle measured counterclockwise in the x-y plane starting from the x axis 

k — wave number 

w = angular frequency 

R = distance from the geometrical center to a certain point in space 
H = directional factor 

F = beam pattern = 201og l0 |H| or 101og 1Q H 2 expressed in dB 

All distances are measured from the origin of the coordinate system to the sources 
along the x and y axes. Phase differences are based upon the phase of the first source, 
which is assumed to be zero. 

The goal of this thesis is to develop a convenient methodology for the 
investigation of the functional dependency of the directional factor on the source 
amplitudes (A, B and C), distances from the origin (D lt D 2 , D 3 and D 4 ), phase 
differences (tp 2 , (p 3 and (p 4 ), and 0 and <l> as demonstrated by the five programs: P3D, 
2DICIR, 2DPCIR, 2DIFIX, 2DPFIX. 
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II. THEORY 



A spatial configuration of simple sources, either discrete or distributed, each with 
its own complex source strength, can be used to represent the quadruplet source case. 
The pressure at a field point is the sum of pressures produced by the individual sources 
[Ref 7: p. 169]. Before we develop the quadruplet source case, the acoustic doublet 
case needs to be described to provide a background for the more complicated 
quadruplet case. 

A. DOUBLET SOURCES 

A doublet source consists of two point sources of strengths of 1 and A, separated 
from the origin, by distances Dj and D 2 and vibrating at the same frequency but with 
a phase difference (p. 




Figure 2.1 Geometry of doublet. 
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( 2 . 1 ) 

(2.2) 



From the geometry shown in Figure 2.1, we have 

cos0 = Sj/Dj, 
sin0 = Ej/Dj, 

rj = V(R * £j ) 2 + V 

= V(R - DjSin©) 2 + (DjCos©) 2 , 

= V(1 - Djsin0/R) 2 + (6jCos0/R) 2 . 

R 

sin0 = c 2 /D 2 , 
cos0 = 5 2 /D 2 , 

r 2 = V(R + 2 + ^ 2 2 

= V(R + D 2 sin0) 2 + (D 2 cos0) 2 , 

— 2- — V(1 + D.,sin0/R) 2 + (5 2 cos0/R) 2 . 

R 

The pressure at point (R,0) due to source l is 
P , - — eK wt ' >“1) 

r l 

where w is the angular frequency. That due to source 2 is 



(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 

(2.3) 



(2.9) 



P = A e j(wt - kr 2 - <p) ( 2 .10) 

r 2 

The total acoustic pressure at point (R,0) is then the sum of P t and P 2 , 

P = JL e i('vt • kR) { JL e ]k(R - r,) + AR e .j[k(r 2 - R) + <p| } (2 1 1) 

R r, r 2 

where R is the distance from the field point to the origin. From (2.4) and (2.S), we get 
k(R-rj) = kR(l - r A /R) 
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( 2 . 12 ) 



= kR[l - V(1 - DjSinG/R) 2 + (Djcos0/R) 2 ], 
k(r 2 - R) = kR(r 2 /R - 1) 



= kR[V(D,sin0/R + l) 2 + (D,cos0/R) 2 - 1 ]. 



(2.13) 



Therefore, equation (2.11) can be rewritten as 



p = _L e j(wt - kR) { JL e jkR[l - V(1 - DjSinG/R) 2 + (DjCosG/R) 2 ] 



R 



+ J^ e -j[kR(V(D 2 sin0/R + l) 2 + (D 2 cos0/R) 2 - 1 ] + <p] } (2 



14) 



In the most frequently encountered cases, the observation of the pressure is made 
at distances greater than the separation of the sources. Therefore, it will be useful to 
derive the far field (R> > (Dj + D-,)) pressure equation from the above equation (2.14). 
From the geometry shown in Figure 2.2, we have 



sin0 = Cj/Dj = £ 2 ,'D 2 . 



(2.15) 



Therefore, (R - rj) and (r 2 - R) can be replaced with DjSinG and D 2 sin0. Also we have 



R R 



(2.16) 



so that the pressure equation becomes 



p _ J_ e j(wt - kR)^ e jkD 1 sin0 + Ae -j(kD 2 sin0 + <p)j 



(2.17) 



|P| = P(R)H(0) 
where P(R) = (1 + A)/R. 



(2.18) 
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H(0) = 



1 



(1 + A) 



^(kDjSinG) + ^ e -j(kD 2 sin0 + <p) 



(2.19) 



This directivity factor, H(0) can be written as 



H(0) = 



1 



(1 + A) 



•{[cos(kD 1 sin0) + Acos(kD 2 sin0 + (p)] 2 



+ [sin(kDjSin0) - Asin(kD 2 sin0 + <p)r) 



2 > 14 



( 2 . 20 ) 




Figure 2.2 Geometry of doublet (far field). 
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B. QUADRUPLET SOURCES 

Similarly, a quadruplet source case can be developed with four sources of 
strengths 1, A, B and C, separated from the origin, by distances Dj, D.,, D 3 and D 4 
and vibrating at the same frequency with phase differences <p 2 , cp 3 and <p 4 (Figure 2.3). 




Figure 2.3 Geometry of quadruplet. 
From the geometry shown in Figure 2.3, we have 



x = RsinGcosO, 


(2.21) 


y = RsinOsin<J>, 


(2.22) 


z = RcosG, 


(2.23) 


R = Vx 2 + y 2 +z 2 


(2.24) 
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where x, y, z are rectangular coordinates and R is the distance from the field point to 
the origin. 

And then 



rj = V(y + Dj) 2 + x 2 + z 2 . (2.25) 

By Taylor's approximation 

r : ~ R + Djy/R + Dj 2 /2R. (2.26) 

rj - R = DjV/R = DjSinGsin^, (2.27) 

r 2 = V(y - D 2 ) 2 + x 2 + z 2 - R - D 2 y/R, (2.28) 

R - 2 ~ D 2 y/R = D 2 sin0sin<l> (2.29) 



where R is much greater than Dp and D 2 - 
Using the same approach 



r 3 = V(x - D 3 ) 2 + y 2 + z 2 = VR 2 - 2D 3 x + D 3 2 

= RVl - (2D 3 .x/R 2 - D 3 2 /R 2 ) . (2.30) 

r 3 ~ R[1 - (1/2)(2D 3 x/R 2 - D 3 2 /R 2 )] = R - D 3 x/R + D 3 2 /2R, (2.31) 

R - r 3 = D 3 x/R = D 3 sin0cos<D (2.32) 

r 4 = V(x + D 4 ) 2 t y 2 4- z 2 ~ R + D 4 x/R, (2.33) 

r 4 - R = D 4 sin0cosO (2.34) 



where R is much greater than D 3 , and D 4 . 

Therefore, the pressure at (R, 0, (p) is 

p = _L e j(wt - kR - kDjSinGsin^) + ^ e j(wt - kR + kD 2 sin0sin<I> - <p,) 

R R 

+ £ e j(wt - kR + kD 3 sin0cos<I> - (p 3 ) 

R 

+ £ e j(wt - kR - kD 4 sin0cosO - (p 4 ) 

R 



(2.35) 
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The directivity factor, 



H(0,O) 



i I e -j(kDjSin0sinO) + ^ e j(kD 2 sin0sinO - (p 2 ) 

(1 + A+B4-C) 

+ g e j(kD 3 sin0cosO - (p 3 ) + £ e -j(kD 4 sin0cosO + cp^ | 



(2.36) 



Expanding the complex terms yields 



H(0,O) 



{[coslkD^inBsinO) + Acos(kD 2 sin0sinO - (p 2 ) 

+ Bcos(kD 3 sin0cos<I> - <p 3 ) + Ccos(kD 4 sin0cosO + <p 4 )] 2 
+ [sin(kDjSin0sinO) - Asin(kD 2 sin0sinO - <p 2 ) 

- Bsin(kD 3 sin0cosO - (p 3 ) + Csin(kD 4 sin0cosO 4- <p 4 )] 2 }' /2 . 



(2.37) 



This is the expression used for the computer program. 

C. GRAPHICAL PRESENTATION OF THE DIRECTIVITY 

1. Two-Dimensional Rectangular and Polar Plots 

From equation (2.37), H (0,0) versus 0, and H (0,0) versus O can be plotted 
if either O or 0 respectively- is given constant value. In a rectangular plot, the 
directivity, F, which is equal to 201og i0 |H (0,0) | is plotted on one axis and 0 or O on 
the other. An example of this is seen in Figure 3.3. In a polar plot, the directional 
factor or beam pattern is a magnitude of the radial distance from the center of the plot 
while the appropriate angle is measured clockwise around the origin. An example of 
this is seen in Figure 3.4. 

2. Three-Dimensional Cartesian Plot 

From equation (2.37), H can be calculated and plotted as a function of both 0 
and O. This provides the three-dimensional surface with its elevation equal to 
201og 10 |H (0,O)| for each 0-0 pair. An example of this is Figure 3.7. This form of 
presentation is difficult to interpret but can be a very useful analytical tool. Further 
explanations are given in the next chapter. 

3. Contour Plot 

A contour plot gives lines of constant H (0,0). As an example, let kDj = 
kD 2 = kD 3 = kD 4 = kD, cp 2 = 0, (p 3 = <p 4 = n and A = B = C = 1. Equation 
(2.37) then can be simplified to 



17 



H(0,O) = 



2 



|cos(kDsin0sin<I>) - cos(kDsin0cos<I>)|. 



(2.38) 



By setting H(0,<I>) = Q, the desired value, and sin0sin<t> = Cj, and sin0cos<D = c 2 
where JCj [ * [c 2 |, we get 



This transcendental equation is then solved for Cj and c,. Now, with Cj and c 2 
determined, we need to solve the second transcendental equation 



which was derived by dividing Cj by c„ for <I> and then, for 0 using the 0-(J> 
relationship. However, many of Cj-c., pairs are possible solutions to the first 
transcendental equation. Therefore, equation (2.40) must be solved for all the different 
Cj-c, pairs by repeating the above procedure. This would result in a contour plot such 
as Figure 3.10. This analytical approach was not pursued further in this research 
because other algorithms for extracting constant values of F for varying 0 and had 
already been obtained and successfully applied. However, the utility of the above 
transcendental equation approach may well bear further investigation, and we will 
pursue a little case-by-case analysis for some cases for which the solution of equations 
(2.39) and (2.40) is relatively straightforward. 



|cos(kDcj) - cos(kDc,)| = 2Q. 



(2.39) 




(2.40) 
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III. PROGRAM DEVELOPMENT 



In analyzing the directional and frequency properties of shaded and phased 
simple arrays several types of display methods were needed. At the first stage of 
developing the computer program, a HEWLETT-PACKARD 86B computer and a 
HP709QA MEASUREMENT PLOTTING SYSTEM were used. By using this 
equipment we were able to get two-dimensional plots as shown in Figures 3.1 and 3.2. 
However, it took much time in executing the program and also we could not do three- 
dimensional plottings. 




Figure 3.1 2-D rectangular plot (HP 86B). 
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PHI VS DIRECTIVITY 



180 (DEG) 




Figure 3.2 2-D circular plot (HP 86B). 

To get a complete understanding of the beam pattern, three-dimensional displays 
viewed from various angles were essential. Consequently, alternative methods of 
display were programmed in FORTRAN and executed on the IBM 3033. There are 
several reasons why this proved to be the most desirable method. Primarily, a three- 
dimensional graphics software package (DISSPLA) is available for use on the IBM 
3033. Secondly, FORTRAN generated the numerical output with greatly increased 
efficiency when using the EXTENDED FORTRAN COMPILER [Ref. 8: p.18]. 
Finally, the IBM 3033 has many more output devices such as VERSTEC, TEK618 and 
SHERPA to name a few. Therefore, the programs used in the HEWLETT- 
PACKARD 86B were rewritten in FORTRAN. The translation from BASIC to 
FORTRAN required some changes in variable names along with the use of some 
convenient function calls. For reference, some of the original BASIC programs are 
attached in Appendix A. 
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Figure 3.3 2-D rectangular plot (DISSPLA). 

In the trial of analyzing only with numerical outputs, BASIC had become 
ineffective in recognizing details due to the very large number of computed values. At 
this stage the graphic package DISSPLA, available on the IBM 3033 computer, was 
used in conjunction with these programs. The DISSPLA package is a library of 
FORTRAN subroutines that facilitate data plotting. It is also device independent, 
w’hich means that it does not rely upon features particular to any type of graphical 
device and has its own extensive symbol and character generation routines. By being 
device independent, programs of the HEWLETT-PACKARD 86B could be translated 
without difficulty. Therefore, DISSPLA was used to produce both two-dimensional 
and three-dimensional graphs. The data used to compute graphs could be generated 
within the FORTRAN program or read from a data file. The former method was used 
here because data files took a great amount of file space due to the large volume of 
values. 
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PHI VS. DIRECTIVITY 



0 



THETA=35 DEGREE 
K=10.5 

Dl=0.5 D2=0.5 
D3=0.5 D4=0.5 
P2= 0 P3= 0 P4= 0 
A=l. B=l. C=l. 




Figure 3.4 2-D circular plot (DISSPLA). 

Two-dimensional graphs were initially tried. By holding one parameter constant 
and with the use of some array variables, plots of directivity versus angle were 
obtained. Examples of two kinds of displays are shown in Figures 3.3 and 3.4. 
Programs are attached in Appendix B. 

To predict the real shape of a directivity pattern in three-dimensions, three- 
dimensional graphics had to be developed. In trying to get some satisfying three- 
dimensional plotting, all methods in DISSPLA were tried. They are SURFUN, and 
SURMAT for regular matrix, scattered points, and three-dimensional vector drawing. 

SURFUN is a subroutine that draws the surface from points in the z direction in 
conjunction with a grid in the x-y plane on the condition that z=f{x,y) is given. One 
constraint that was faced in applying this function was that there was more than one z 
value for a specific x-y coordinate. 
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Figure 3.5 3-D plot(SURMAT for scattered points). 

The next trial was SURMAT for a regular matrix (for example, a rectangular 
three-dimensional lattice of points). For the same reason as in applying the SURFUN 
subroutine, this also turned out to be non-applicable even though all data points were 
translated into a regular matrix. 

The next option tried was SURMAT for scattered points (for example, a three- 
dimensional irregular cluster of points). As seen in Figure 3.5 it also was unsuccessful 
due to more than one point on one specific coordinate of the x-y plane. The program 
is listed in Appendix C for reference. 

The last option, the three-dimensional vector drawing method gave us a more 
understandable display as shown in Figure 3.6. However, it did not show the beam 
pattern clearly because all the lines were visible, even from different view angles. At 
this stage the problem was the removal of hidden lines. Unfortunately, that was not 
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possible because the starting and ending points of the hidden lines could not found. 
The program is listed in Appendix D. 




Figure 3.6 3-D plot(vector drawing method). 

Consequently, with the consideration of all the constraints mentioned above, 
another type of format was made to see the details of the beam pattern. Now, a plot 
was made using the function SURMAT about the 0-0 plane instead of the x-y plane. 
This eliminated the previous multiple value problem. The results showed the most 
desirable presentation among all those tried. As shown in Figure 3.7 it gives a clearer 
appreciation of the directivity pattern in conjunction with 0 and O. The program is 
listed in Appendix E. 

We can relate this three-dimensional graph (Figure 3.7) to the two-dimensional 
graphs (Figures 3.3 and 3.4) for interpretation. As shown in Figure 3.3 the x-axis is 
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the angle and the y-axis is the directivity which is expressed in dB. Figure 3.4 is the 
circular graph with angle working counterclockwise in 10° increments and each circle is 
an increment of 10 dB whth the center of - 40 dB. Figure 3.7 is the three-dimensional 
graph from the view point of 20° from the x-z plane, and 60° from the x-y plane. Here, 
the x-axis is <I>, the y-axis is 0 and the z-axis is the beam pattern (dB). From Figure 
3.3 four major points can be found at approximately 50°, 130°, 230° and 310° of 0 
with O fixed at 126°. From Figure 3.7 the same feature can be determined at 50° of 0, 
126° of 0 as seen in Figure 3.3. Since Figure 3.3 was made by holding <I> constant at 
126°, one can look along the <I> = 126° line in Figure 3.7 and note the null at 0 = 
90°, the maximum at 0 = 50°, and the relative maximum at 0 = 0°. Note the heavy 
lines in Figure 3.7. 



DIRECTIVITY PATTERN 




K=10.5 

Dl= 0.5 D2= 0.5 
D3= 0.5 D4= 0.5 
VIEW ANGLE=+20,+60 
P2=0 P3= 0 P4= 0 
A=L B=i. C=l. 



Figure 3.7 3-D plot(SURMAT in 0-<D plane). 
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From the circular graph (Figure 3.4) four major points can be found approximatly at 
40°, 140°, 220° and 320° of <I> with 0 fixed at 35°. In the three-dimensional graph one 
of these major points is seen approximately at 35° of 0, 40° of <J> and the other major 
point is seen approximately at 35° of 0, 140° of <I>. Comparison is made easier when 
we see Figure 3.8 and Figure 3.9 which were drawn with the surfaces cut at O = 126° 
and 0 = 35° respectively. Note the shift in axes between the two figures. All these 
points of two-dimensional graphs correspond to those of the three-dimensional graph. 




Figure 3.8 3-D plot(SURMAT in 0-<P plane). 

Inversely, every point of the three-dimensional graph can be read in two-dimensional 
graphs to see detailed values whenever needed. 

Finally, a trial was made to draw a constant value plot of the beam pattern. The 
first approach was to derive the functions of <J> = f (0) or 0 = f (<J>) and plot these on 
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DIRECTIVITY PATTERN 




K=10.5 

Dl= 0.5 D2= 0.5 
D3= 0.5 D4= 0.5 
VIEW ANGLE=+20,+60 
P2= 0 P3= 0 R=0 
A=l. B=l. C=l. 



Figure 3.9 3-D plot(SURMAT in 0-<P plane). 

the 0-<I> plane. However, failing to derive these relationships for the general four 
element case made this method of making a contour plot of a fixed F value impossible. 
Instead, a modified version of the existing 0-<I> program was used. It sorted on a small 
range of values about the desired value of F to find the corresponding 0-<l> coordinates 
which were then plotted as in Figure 3.10. This type of plot is included throughout the 
study. 

With these forms of display and the speed of the FORTRAN version of this 
program we can now use this method to study the directional and frequency properties 
of shaded and phased simple arrays. 
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IV. COMPARISON WITH SAMPLE CASES 



In making comparisons, some simple parameters were used to make the 
comparisons easier. There are three cases included in this chapter. Each case has a 
unique set of input data, and computer generated and experimental plots. In all cases, 
a two points were picked and hand calculations made to verify the computer output. 
An example of the hand calculations is given for the first case only. Comparisons were 
also made between hand or computer generated values and the two-dimensional polar 
plot, between the theoretical two-dimensional plot and the experimental data plot, 
between the tw'o-dimensional plot and the three-dimensional plot, and finally between a 
fixed value of the beam pattern (a "track") and its three-dimensional plot. 

While in general, 0(0) or 0(0) for a constant F could not be obtained 
analytically, in a few cases, included below, it was possible to obtain the functional 
dependencies. These are discussed in the following cases. 

A. CASE 1 

When kD T = kD-> = kD 3 = kD 4 = 5.275, 

<p 2 = 0, q> 3 = (p 4 = rt, 

A = B = C = 1 

equation (2.37) can be simplified as 



H(9,<D) 



|cos(5.275 sin0sinO) - cos(5.275 sin0cosO)|. 



(4.1) 



Using the above input data with values of 0 = 40° and <J> = 0° we get H(0,O) = 
0.984566 and F (0,0) = -0.135104 from the computer output as seen in Appendix G. 
If we use the above simplified equation and hand-computation, we get H(0,O) = 
0.98457 and F(0,O) = -0.135068 . 

In the same way, if 0 = 90° and <J> = 0° then we get H = 0.233302, and F = 
-12.641622 from the computer output and H = 0.23330, and F = -12.6417 from the 
hand-calculation. These hand-calculated values are almost the same the the computer 
output. Comparisons can also be made between these selected values and the two- 
dimensional plot, Figure 4.1a. Additionally, the above two-dimensional plot 
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corresponds well in shape with Figure 4.1b which is from LCDR John Butler's 
experiment [Ref. 10]. 

If we compare the shape of the two-dimensional plot (Figure 4.1a) with that of 
the three-dimensional plot, Figure 4.1c, we see that beam pattern of the two- 
dimensional graph is identical in the <I> = 0° plane of the three-dimensional plot. 




Figure 4.1a 2-D plot from the computer output. 

A two-dimensional plot of equation (4.1), of F values between -10.9 dB and -10.0 
dB results in Figure 4. Id. For F values less than -90.0 dB we get Figure 4.1e. 

In the case where F < -90 dB, H < 10‘^' ^ which is small enough to be assumed to 
be 0, |c 1 1 ~ |c 2 | or 0 = ntt n = 0,1,2 — , and equation (2.40) can be tanO = 1 or -1. 
From these resulting equations we can extract <I> = 45° + ntt/2, n = 1,2 — . As seen 
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in Figure 4.1e, when F < -90 dB the asymptote matches the nulls in the 0-<P 
relationship of Figure 4.1c. Note the nulls at 0 = 0°, 180° for all <l>, and <l> = 45°. 
These figures of constant dB level aid in the interpretation of Figure 4.1c. The 
modified P3D program (TRACE) is included in Appendix F. 



QUADRUPOLE BEAM PATTERN (dB) 




R-3.20 M kcMO. 55 FREQ-960 Hz 
PHASE-0. 180. 0. 180 Deg AHPL- 0. 265 Vac ANGLE-0 Deg 



Figure 4.1b 2-D plot from the experiment. 
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DIRECTIVITY PATTERN 




K=10.55 

Dl= 0.5 D2= 0.5 
D3= 0.5 D4= 0.5 
VIEW ANGLE=-!-20,+60 
P2= 0 P3=180 P4=180 
A=l. B=l. C=l. 



Figure 4.1c 3-D plot from the computer output. 
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B. CASE 2 

In this case, a triplet array is investigated. Referring to Figure 2.3, the fourth 
element is turned ofF and its opposing element is moved toward the center of the array 
so its distances is half that of the two elements. Also, the source 2 is 180° out of phase 
with the other two elements. 

When kDj = kD 2 = 5.2758, 

kD 3 = 2.6379, kD 4 = 0, 

(p 2 = 7t, q> 3 = (p 4 = 0, 

A = B = 1, C = 0 
the equation (2.37) can be simplified as 

H(0,4>) = — {[cos(2.6379 sinGcosO)] 2 + [2 sin(5.2758 sinGsinO) 

- sin(2.6379 sinGcosO)] 2 } >/z. (4.2) 

If values are calculated from the equation (4.2) as in case 1, then comparisons can be 
made between these calculated values and the two-dimensional plot as seen in Figure 
4.2a (for 0 = 0°, 20° and 90°, and <P = 90°, we get F = -9.5424, -2.7431 and -3.6776 
respectively from hand-computation.). Additionally, the two-dimensional plot 
corresponds well in shape with Figure 4.2b from the experiment (Ignore the ordering of 
the relative phases in this figures; it is the consequence of a different labeling 
convention). 

If we compare the shape of a two-dimensional plot (Figure 4.2a) with that of the 
three-dimensional plot (Figure 4.2c), we see that beam pattern of the two-dimensional 
graph is identical in the <I> = 90° plane of the three-dimensional plot. 

For a contour plot of this case, a two-dimensioal plot of equation (4.2) of F 
values between -5.9 dB and -5.0 dB results Figure 4. 2d. For F values between -10.9 dB 
and -10.0 dB, we get Figure 4.2e. For F values less than -35.0 dB, we get Figure 4.2f. 
From equation (4.2), we get F = -35.876 for 0 = 143° and <I> = 10° from hand- 
calculation. This point is seen in Figure 4.2f. Successive examination of Figure 4.2d 
through 4.2f shows how the contour lines migrate as lower and lower values of F are 
examined. These figures of constant dB level aid in the interpretation of Figure 4.2c. 
The modified P3D program (TRACE) is included in Appendix F. 
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Figure 4.2a 2-D plot from the computer output. 
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TRIPOLE BEAM PATTERN (dB) 




R-3.20 M FREQ-960 Hz 

PHASE-O. 0. 180 Dog AMPL- 1. 066 Vac ANGLE-90 Dog 



Figure 4.2b 2-D plot from the experiment. 
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DIRECTIVITY PATTERN 




K=17586 
Dl= 0.3 D2= 0.3 
D3= 015 D4= 0. 

VIEW ANGL3=+20,-h60 
P2=180 P3= 0 P4= 0 
A=L 3=1. C=0. 



Figure 4.2c 3-D plot from the computer output. 
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TRACK (-5.0 DB) 




Figure 4. 2d Track. (-5 dB). 
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TRACK (-35.0 DB) 




Z = -35.0 DB 
K=17.586 
Dl=0.3 D2=0.3 
D3=0.15 D4=0. 
P2=180 P3= 0 P4= 0 
A=l. B=l. C=0. 



Figure 4.2f Track (-35 dB). 
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C. CASE 3 

In this case, a doublet array is investigated. Referring to Figure 4.2a, this is the 
same as the first case except that the third and fourth elements are turned off and all 
elements are in phase. 

For kD 2 = kD ; = kD 3 = kD 4 = 5.275, 

<p 2 = <p 3 = <p 4 = o, 

A = 1, B = C = 0 
equation (2.37) can be simplified as 

H(0,O) = |cos(5.275 sin6sin<P)| 

As in cases l and 2, for 0 = 0°, 35°, 90° and <I> = 90°, we get F = 0, -0.058 and 
-5.459 respectively from hand-computation. These points are found in Figure 4.3a and, 
the two-dimensional plot again corresponds well in shape with Figure 4.3 b from the 
experiment. 

If we compare the shape of a two-dimensional plot (Figure 4.3a) with that of the 
three-dimensional plot (Figure 4.3c below), we see that the beam pattern of the two- 
dimensional graph is identical in the <I> = 0° plane of the three-dimensional plot. 

For a constant value of H(8,<J>), we can follow the same procedure to get the 
equation for contour plots as in case I and 2. The resultant figures are Figures 4.3d, 
4.3e and 4.3f. Figure 4.3d is drawn for the F values between -5.9 dB and -5.0 dB and 
Figure 4.3e is for the F values between -8.9 dB and -8.0 dB. Figure 4.3f is for F values 
less than -35.0 dB. 
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Figure 4.3a 2-D plot from the computer output. 
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DIPOLE BEAM PATTERN (dB) 




R-3.20 H kd-10. 55 FREQ-960 Hz 
PHASE-0.0 Dog AMPL-6. 53 Vac ANGLE-90 Oog 



Figure 4.3b 2-D plot from the experiment. 
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Figure 4.3c 3-D plot from the computer output. 
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TRACK (-8 DB) 



O 




Z = -8.0 DB 
K=10.55 
Dl=0.5 D2=0.5 
D3=0.5 D4=0.5 
P2= 0 P3= 0 P4=0 
A=l. B=0. C=0. 



Figure 4.3e Track (-8 dB). 
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TRACK (-35.0 DB) 



o 




Z = -35.0 DB 
K=10.55 
Dl=0.5 D2=0.5 
D3=0.5 D4=0.5 
P2= 0 P3= 0 P4= 0 
A=L B=0. C=0. 



Figure 4.3f Track (-35 dB). 



V. EVALUATION 



In the last chapter, we looked at three simple cases. These established the 
accuracy of the five computer programs. For example, in comparisons between hand 
calculation and computer output, both sets of values were proven nearly identical. 
Again, due to the complexity of calculations, not all values were hand-calculated. 

When a two-dimensional circular graph was compared with the selected hand- 
computed values, both groups of values were again virtually identical. As another test 
of the validity of the five computer programs, some of the computer generated two- 
dimensional plots were compared against those from an experiment [Ref. 9: p. 1 8]. 
Visual comparisons of the two types of plots were similar in shape, but values in dB 
were somewhat different. These differences might be caused by the normalization of 
the angular dependance of the sound pressure. Also, because of the directivity of the 
real sources, there is a slight decrease in dB level of the experimental plots when the 
array pointed away from the receiver. However, they are small enough that the 
comparisons still seem quite reasonable. 

Next, two-dimensional versus three-dimensional graphs were compared. 
Agreement here shows the programs are working properly. By varying the range of 0 
and <t> that are plotted, and also by varing the view angles, values and the shape of the 
curve can be visualized without great difficulty. 

Tracks of certain directivity values were drawn using the modified P3D program 
(TRACE) and were compared to the three-dimensional plots. The shape of the 
resulting tracks corresponded well with those of the three-dimensional patterns, again 
showing proper functioning of the programs. 

With all the comparisons made in Chapter IV, the five computer programs 
developed in this thesis were proven effective in predicting the directivity pattern in all 
directions of space with or without differences in the strength of the sources, inter- 
source distances and phases. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The following conclusions are offered. 

• The P2DICIR, P2DPCIR, P2DIFIX, P2DPFIX programs successfully 
predict the directivity beam pattern in two dimensions. 

• The P3D program predicts and visualizes the directivity beam pattern in the 
three-dimension space. 

• They allow the presentation of complete beam patterns, which contain 
considerable information, in a succinct, easy to comprehend and interpret form. 

B. RECOMMENDATIONS 

The following recommendations and areas of further investigation is suggested. 

• Further studies for computer programs to visualize the beam pattern in the 
x-y-z space are needed to eliminate hidden lines in three-dimensional plotting. 

• Transcendental equations need to be solved to get the contour plot for fixed 
F values. 
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APPENDIX A 

BASIC PROGRAMS (HP COMPUTER) 



xxxx x x xx xx xxxxx xxxxxx x * INFORMATION 

X 

* THE OBJECTIVE OF THESE PROGRAMS IS TO CALCULATE AND 

* DRAW THE DIRECTIVITY PATTERN IN TWO-DIMENSION. 

* THIS IS WRITTEN IN BASIC LANGUAGE AND PLOTTER IS HP7090A 

* 

XXXXXXXXXXXXXXXXXXXXXX VARIABLE DEFINITION X X XX XX XXX XXX XXXX XX X* 
* 

* Z : DIRECTIVITY IN dB 

* P2 : PHASE (DEGREE) OF SOURCE A 

* P3 : PHASE (DEGREE) OF SOURCE B 

* P4 : PHASE (DEGREE) OF SOURCE C 

* P5 : VARIABLE TO CHANGE DEGREE TO RADIAN 

* P6 : RADIAL VALUE OF PI 

* P7 : RADIAL VALUE OF P2 

* P8 : RADIAL VALUE OF P3 

* K : WAVE NUMBER 

* D : INTER-SOURCE DISTANCE 

* P : PHI IN DEGREE 

* I : THETA IN DEGREE 

* PI : 3 . 19-15917 
x 



-X 



X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 



xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 



xx x xxxxxxxxxxxx xxxxx xxxxxxxx x xxxxxxx x xx x x xxxxx xxxx x x xx xxxxxxxxx xx x xxxxxx 

HP PROGRAM FOR THE RECTANGULAR GRAPH 

X XXXXX XX XXX X XXXXXXXXX XXXX XX X X XX XX X * X XX XXXXXXXXXXXX X XXXXXXXXXXXXXXXXXX XXX 



INPUT PARAMETERS 



5 OPTION BASE 1 

6 RAD 

10 DIM Z( 360 ) 



20 


DISP "WHAT 


IS 


THE 


VALUE 


OF 


KD? 


n 

o 


KD=K*D 






30 


INPUT KD 




















60 


DISP "WHAT 


IS 


THE 


VALUE 


OF 


P2( DEG )? P2= "! 


P2=PHASE 


OF 


A 


61 


INPUT P2 




















62 


DISP "WHAT 


IS 


THE 


VALUE 


OF 


P3( DEG )? P3 = M ! 


P3=PHASE 


OF 


B 


63 


INPUT P3 




















64 


DISP "WHAT 


IS 


THE 


VALUE 


OF 


P3( DEG )? P4= " ! 


P4=PHASE 


OF 


C 


65 


INPUT P4 




















66 


DISP "WHAT 


IS 


THE 


VALUE 


OF 


A? 


A= " 








67 


INPUT A 




















68 


DISP "WHAT 


IS 


THE 


VALUE 


OF 


B? 


B= " 








69 


INPUT B 




















70 


DISP "WHAT 


IS 


THE 


VALUE 


OF 


c? 


C= " 








71 


INPUT C 




















72 


DISP "IS P 


FIXED? 


( YES jAN = : 


L > 


NO jAN=0 ) AN= " 









73 INPUT AN 

74 IF AN=1 THEN DISP "WHAT IS THE VALUE OF P( DEG )? P= "ELSE GOTO 76 

75 INPUT P 

76 DISP "IS I FIXED? (YESjANS=l > NO>ANS=0) ANS= " 

77 INPUT ANS 

78 IF ANS=I THEN DISP "WHAT IS THE VALUE OF I (DEG)? 1= "ELSE GOTO 82 

79 INPUT I 

83 IF AN=0 AND ANS=0 THEN GOTO 86 

84 IF AN=0 AND ANS=1 THEN GOTO 90 

85 IF AN=1 AND ANS=0 THEN GOTO 93 

86 DISP "x" *"KxD=" *KD >" 

87 DISP "*PHASE OF A( DEG )=" *P2)"PHASE OF B< DEG )=" >P3 >"PHASE OF C( DEG )=" JP4>"* 

88 DISP "*">"A="*A>" B=" JB > "C=" jC > £ GOTO 99 

90 DISP "*"j"K*D="‘>KD >" I(DEG)="jI>" ">"x" 

91 DISP "*PHASE OF A( DEG )=" jP 2* "PHASE OF B( DEG ) = " jP3 j "PHASE OF C( DEG )=" jP4 

92 DISP j"A=" ; A>" B= n 5B,"C=">C,"*" 3 GOTO 103 

93 DISP "*'V , K^D= ,, jKD," I( DEG )=" jl » " ","*" 
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********* 



94 OISP "*PHASE OF A( DEG )=">P2 $ "PHASE OF B( DEG )='SP3 $ "PHASE OF C( DEG )='SP4,"*" 

95 DISP 

100 DISP " " 

101 DISP " ">" THETA'S" PHI"," DIRECTIVITY" 

102 FOR P = 1 TO 360 0 GOTO 110 

103 DISP “ XXXXXXXXXXXXXXXXXXXX XX X X XXXXXXXXXXXXXXXXXXX XXX X XXXX XX " 

104 DISP " " 

105 DISP " ">" THETA'S" PHI'S" DIRECTIVITY" 

106 FOR P=1 TO 360 0 GOTO 115 

107 DISP " XXX -X XXXXXX XX XXXXXXXX XX X X XXXXXXXXXXXXXXXXXXXXXXXXXXXXX 11 

108 DISP " " 

109 DISP " THETA'S" PHI'S" DIRECTIVITY" 



MAIN PROGRAM 



110 FOR 1=1 TO 360 

115 P5=PI /180 

116 P6=P2*P5 3 P7=P3*P5 3 P8=P4xPS 

117 Q=P-90 

118 Q1=SIN ( Q*P5 ) 

119 I1=SIN ( I*P5 ) 

120 Q2=SIN ( P*P5 ) 

121 Wl=KD*Il*Q2/2 
121 W2=KD*Il*Ql/2 

127 W=COS (W1)+A*C0S < W1+P6 )+B*COS ( W2-P7 )+C*COS (W2+P8) 

128 G=SIN (W1)-A*SIN ( W1+P6 )+B*SlN ( W2-P7 )-C*SlN (W2+P8) 
130 H=SQR ( W**2+G**2 )/( 1+A+B+C ) 

140 Z(P)=20*LGT <H) 

143 T=20XLGT (ABS (SIN t.5*KD*Iim 
150 DISP USING 153 j I»P,Z(P) 

153 IMAGE 5X»3D . D >5X >3D . D > 15X >SDDD . DDD > 15X,SDDD . DDD 

155 ! IF Z( IK-20 THEN I=I+.l 0 GOTO 119 

156 ! 1= IP(I) 

157 ! IF Z( P K-20 THEN P=P+.l 2 GOTO 117 

158 ! P= IP( P)) 

160 IF AN=1 AND ANS=0 THEN GOTO 163 ELSE GOTO 161 

161 IF AN=0 AND ANS=1 THEN GOTO 164 ELSE GOTO 163 

162 NEXT P 0 GOTO 96 

163 > NEXT I 

164 NEXT P 



PLOTTING WITH HP7090A 



170 PLOTTER IS 705 
180 CLEAR 
190 GCLEAR 
200 GRAME 

210 LOCATE 15,145,10,80 

220 SCALE 0,390,-240,0 

230 CLIP 0,360,0,-240 

240 AXES 5,5,0,-240,10,10,2 

250 LINE TYPE 1 

260 GRID 20,-20 

270 FRAME 

280 LORG 6 

290 CSIZE 4 

300 FOR 1=0 TO 360 STEP 40 
310 MOVE 1,-241 
320 LABEL I 

330 NEXT I 
340 LORG 8 

350 FOR I=-240 TO 0 STEP 40 

360 MOVE -1,1 

370 LABEL I 

380 NEXT I 

390 CSIZE 5 

400 LORG 5 

410 MOVE 180,40 

420 CSIZE x 

430 PEN 2 

440 LABEL "PHI VS DIRECTIVITY" 

450 LORG 5 

460 MOVE 180,20 

470 CSIZE 4 

480 LABEL "(THETA = 45 DEGREE)" 

490 CSIZE 5 

500 LORG 5 

510 MOVE 180,-26- 

520 CSIZE 4 
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530 LABEL "PHI (DEG)" 

540 LOIR 90 
550 MOVE -30,-120 
560 DEG 
570 LDIR 90 
580 CSI2E 4 

590 LABEL "DIRECTI VITY ( dB ) M 

600 CSI2E 5 

610 MOVE 0,0 

620 FOR P=1 TO 360 

630 DRAW P >2( P ) 

640 NEXT P 
650 PEN 0 
700 END 
117392 



HP PROGRAM FOR CIRCULAR GRAPH 



INPUT PARAMETERS 



5 OPTION BASE 1 

6 RAD 

10 DIM X( 360 ) 

11 DIM Y ( 360 ) 

12 DIM 2(360) 



20 


DISP "WHAT 


IS 


THE 


VALUE 


OF 


KD? 


KD= " ! 


KD=K«D 






30 


INPUT KD 




















60 


DISP "WHAT 


IS 


THE 


VALUE 


OF 


P2( DEG )? P2= "! 


P2=PHASE 


OF 


A 


61 


INPUT P2 




















62 


DISP "WHAT 


IS 


THE 


VALUE 


OF 


P3( DEG )? P3= " ! 


P3=PHASE 


OF 


B 


63 


INPUT P3 




















64 


OISP "WHAT 


IS 


THE 


VALUE 


OF 


P4( DEG )? P4= " ! 


P4=PHASE 


OF 


C 


65 


INPUT P4 




















66 


DISP "WHAT 


IS 


THE 


VALUE 


OF 


A? 


A= " 








67 


INPUT A 




















68 


DISP "WHAT 


IS 


THE 


VALUE 


OF 


B? 


B= " 








69 


INPUT B 




















70 


DISP "WHAT 


IS 


THE 


VALUE 


OF 


C? 


C= " 








71 


INPUT C 




















72 


DISP "IS P 


FIXED? 


( YES ;AN=: 


L , 


NO ; AN=0 ) AN= " 









73 INPUT AN 

74 IF AN=1 THEN DISP "WHAT IS THE VALUE OF P( DEG )? P = "ELSE GOTO 76 

75 INPUT P 

76 DISP "IS I FIXED? (YES)ANS=1 , NO;ANS=0) ANS= " 

77 INPUT ANS 

78 IF ANS=1 THEN DISP "WHAT IS THE VALUE OF I( DEG )? 1= "ELSE GOTO 82 

79 INPUT 1 

83 IF AN=0 AND ANS=0 THEN GOTO 86 

84 IF AN=0 ANO ANS=1 THEN GOTO 90 

85 IF AN=1 AND ANS=0 THEN GOTO 93 

86 DISP )"K*0=")KD," 

87 DISP "*PHASE OF A( DEG )=" ) P2 \ "PHASE OF B( DEG )=")P3) "PHASE OF C( DEG )=" ;P4,"*" 

88 DISP 5 "A = " ; A , " B=" jB , "C=" )C , 3 GOTO 99 

90 DISP "*" )"K*0=")KD," I ( DEG )=" jl 

91 OISP "*PHASE OF A( DEG )=" )P2 )" PHASE OF B( DEG )=" )P3 ; "PHASE OF C( DEG )=" )P4>"*" 

92 DISP "*" ; "A=" ) A, " B=" ;B, "C=" ;C >"*" 3 GOTO 103 

93 DISP )"K*0=" jKD ," I(0EG)=")I," 

94 DISP "*PHASE OF A( DEG )=" ;P2 \ "PHASE OF B( DEG )=" \ P3 )"PHASE OF C( DEG )=" )P4, 

95 DISP "*"5"A=";A," B=" 5 B,"C=" jC,"*" 3 GOTO 107 

99 DISP 

100 DISP " " 



101 DISP 



V THETA";" 



102 FOR P=1 TO 360 3 GOTO 110 

103 DISP "**** 

104 DISP " " 



PHI"," 



DIRECTIVITY" 






105 DISP " ";" THETA* V 



PHI* 



DIRECTIVITY" 



106 FOR P=1 TO 360 3 GOTO 115 

107 DISP " XXXXXXXXttXXXX* * **X * KXXXXXXM X* » X XXXX*XXXXXXX * XXX*XXXXX M 

108 DISP " " 

109 DISP " ")" THETA")" PHI"," DIRECTIVITY" 



MAIN PROGRAM 



110 FOR 1=1 TO 360 
115 P5=PI /180 
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116 P6=P2*P5 3 P7=P3*P5 2) P8=P4*P5 

117 Q=P-90 

118 Q1=SIN (Q*P5) 

119 I1=SIN <I*P5) 

120 Q2=SIN ( P*P5 ) 

121 Wl=KD*Il*Q2/2 
121 W2=KD*Il*Ql/2 

127 W=COS IW1HA*C0S t W1+P6 HB*COS ( W2-P7 )+C*COS (W2+P8) 

128 G=SIN (W1)-A*SIN t W1+P6 )+B*SIN ( W2-P7 )-C*SIN (W2+P8) 
130 H=SQR ( W**2+G**2)/l 1+A+B+C) 

140 2( I ) = 20*LGT (H3+240 

141 X(I)=Z(I)*COS ( I*P5 ) 

142 Y( I )=Zt I )*SIN i I*P5 ) 

143 T=20*LGT (ABS i SIN (.5*KD*Il)n 
150 DISP USING 153 > I,P,Z(P) 

153 IMAGE 5X,SDDD.DDD,15X>SDDD.DDD 

1 55 ! IF Z( IK-20 THEN I=I+.l D GOTO 119 

156 ! 1= IP(I) 

157 ! IF Z( P K-20 THEN P=P+.l 2> GOTO 117 

158 ! P= IPCP)) 

160 IF AN=1 AND ANS=0 THEN GOTO 163 ELSE GOTO 161 

161 IF AN=0 AND ANS=1 THEN GOTO 164 ELSE GOTO 163 

162 NEXT P 2) GOTO 96 

163 NEXT I 

164 ! NEXT P 



PLOTTING WITH HP7090A 



170 


PLOTTER IS 705 




180 


CLEAR 




190 


GCLEAR 




200 


GRAME 




210 


LOCATE 25,120,5,80 




220 


SCALE -260,320,-220 


,240 


230 


CLIP -240,240,-240, 


240 


235 


XAXIS 0,20,-240,240 




240 


YAXIS 0,20,-240,240 




250 


LINE TYPE 1 




270 


FRAME 




280 


LORG 6 




290 


CSIZE 4 




300 


FOR 1=0 TO 360 STEP 


40 


310 


MOVE 1,-5 




320 


LABEL 1-240 




321 


MOVE -1,1 




322 


LABEL 1-240 




330 


NEXT I 




340 


LORG 6 




341 


MOVE 0,-5 




342 


LABEL -240 




350 


FOR I=-240 TO -60 STEP 60 


360 


MOVE -1,1 




370 


LABEL -240-1 




371 


MOVE 1,-5 




372 


LABEL -240-1 




375 


NEXT I 




380 


MOVE 0,0 




381 


FOR 1=0 TO 240 STEP 


60 


382 


FOR J=0 TO 360 




383 


X=I*COS (J*PI /180 ) 


384 


Y=I*SIN (J*PI /180 ) 


385 


DRAW X ,Y 




386 


NEXT J 




387 


NEXT I 




390 


CSIZE 5 




400 


LORG 5 




410 


MOVE 0,310 




420 


CSIZE 8 




430 


PEN 2 




440 


LABEL "THETA VS DIRECTIVITY' 


450 


LORG 50 




460 


MOVE 0,285 




470 


CSIZE 4 




480 


LABEL " ( PHI = 90 DEGREE )" 


490 


CSIZE 5 




500 


LORG 5 




510 


MOVE 300,0 




520 


CSIZE 4 




530 


LABEL M 0( DEG )" 




550 


MOVE -300,0 




580 


CSIZE 4 




590 


LABEL "180( DEG ) M 
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610 MOVE 0,0 
620 FOR 1-1 TO 360 
630 DRAW X( I ),Y( I ) 
640 NEXT I 
650 PEN 0 
700 END 
111174 
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APPENDIX B 

P2DIFIX, P2DPFIX, P2DICIR AND P2DPCIR 



x x xxxxx xxxxxx xxx xxxxxx x INFORMATION xxxxxxxxxxxxxxxxxxxxxxxx »»»*»*x»*» 

X X 



x THE OBJECTIVE OF THESE PROGRAMS IS TO CALCULATE AND * 

* DRAW THE DIRECTIVITY PATTERN IN TWO-DIMENSION. * 

x THE. DRAWING METHOD IS THE SUBROUTINE CURVE OF DISSPLA. * 

x 1(0) AND P(<P) CAN BE USED INTERCHANGEABLY FOR * 

* P2DPFIX AND P2DICIR PROGRAMS, 
x 

xxxxx x xx x x x x x x xxx x xxxx VARIABLE DEFINITION 
x x 



X 


Z : 


DIRECTIVITY IN dB 












X 


X 


PI 


: PHASE (DEGREE) OF SOURCE A 












X 


X 


P2 


: PHASE (DEGREE) OF SOURCE B 












X 


X 


P3 


: PHASE (DEGREE) OF SOURCE C 












X 


X 


P4 


: VARIABLE TO CHANGE DEGREE TO RADIAN 










X 


X 


P5 


: RADIAL VALUE OF PI 












X 


X 


P6 


: RADIAL VALUE OF P2 












X 


X 


P7 


: RADIAL VALUE OF P3 












X 


X 


K : 


WAVE NUMBER 












X 


X 


D1 


: DISTANCE FROM THE GEOMETRICAL 


CENTER 


TO 


THE 


SOURCE 


1 


X 


X 


D2 


: DISTANCE FROM THE GEOMETRICAL 


CENTER 


TO 


THE 


SOURCE 


A 


X 


X 


D3 


: DISTANCE FROM THE GEOMETRICAL 


CENTER 


TO 


THE 


SOURCE 


B 


X 


X 


04 


: DISTANCE FROM THE GEOMETRICAL 


CENTER 


TO 


THE 


SOURCE 


C 


X 


X 


P : 


PHI IN DEGREE 












X 


X 


I : 


THETA IN DEGREE 












X 


X 


PH 


: PHI IN RADIAN 












X 


X 


TH 


: THETA IN RADIAN 












X 


X 


PHI 


: DUMMY VARIABLE OF PHI 












X 


X 


THETA : DUMMY VARIABLE OF THETA 












X 


X 


PI 


: 3.1415917 












X 


X 




• 












X 



xxxxxxxxxxxxxxxxxxxxxx VARIABLE DECLARATION **** 

C 

C 

xxxxxxxxxxxxxxxxxxxxx xx x xx x xxxxx xx xxxxxxxxxxxxxx x x 



PROGRAM FOR THE TWO-DIMENSIONAL RECTANGULAR GRAPH 
xxxxxxx xxxx xx xxxx xx x xx xxxxx xxxx xx xx xxxxx xx xxx xx xxxx x x xxxxxx xxxx x xx x xx x xx 



c 

DIMENSION PHI( 361 ) ,THETA( 361 ) ,Z( 360 ) 

REAL P4,P5,P6,P7,Wl,W2,W3>W4,W,G,H,K>Dl,D2>D3,D4>PH,TH>PI>a,b,c 
INTEGER P,P1,P2,P3,I 
C 

C 

C INPUT PARAMETERS 

C 

c 

PI=3. 1415917 
K=5 .5 
Dl=l. 

D2=l. 

D3 = l. 

D4=l. 

Pl= 180 
P2= 0 
P3= 180 
A = 1. 

B = 1. 

C = 1. 

I = 25: USE P FOR P2DPFIX 
C 

C 

C MAIN PROGRAM 

C 

C 

P4=PI/180 

P5=P1*P4 

P6=P2XP4 

P7=P3*P4 

c DO 100 1=1,360: MOVE C TO LINE NO. 76 FOR P2DPFIX 
DO 200 P=l,360 
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* * 



ooono o * * o n n oaooo 



TH=I*P4 

PH=P*P4 

W1=K*D1*SIN( PH )*SIN< TH ) 

W2=K*D2*SIN! PH )*SIN( TH )_P5 
W3=K*D3*SIN! TH )*COSl PH )-P6 
W4=K*D4*SIN! TH )*COSl PH )+P7 
H=COS( W1 )+A*COSl W2 )+B*COS! H3 )+C*COSl W4 ) 
G=SINl W1 )-A*SINl W2 )-B*SIN( W3 )+C*SINl W4 ) 
H=SQRT( W*W+G*G )/( 1+A+B+C ) 

Z( P )=20*ALOG10( H ) 

IF! Z! P ) . LE . -40 . ) THEN 
Z( P )=-40 

ELSE 

Z!P)=Z!P) 

ENO IF 
THETA! P)=I 
PHI ! P )=P 

WRITE ( 8 >1 ) PHI!P),Z!P) 

1 FORMAT! 2X,2I E10.3,5X) ) 

200 CONTINUE 

100 CONTINUE 



A PROGRAM FOR PLOTTING BY SHERPAOR BY TEK618 



CALL MEDBUF 
CALL COMPRS 

CALL TEK618 : MOVE C TO ORAW USING TEK618 
CALL PHYSORd. ,1. ) 

CALL PAGE( 15.5,12. ) 

CALL AREA2D! 10.5,9.1) 

CALL FRAME 
CALL COMP LX 
CALL HEIGHT ( .2) 

CALL XNAME! 1 PHI( DEG) $' ,10) 

CALL YNAME( ' DIRECTIVITY! DB )$ 1 > 15 ) 

CALL YTICKS( 8 ) 

CALL XTICKS! 8 ) 

CALL GRAF( 0 ,40,360,-40 ,4,0 ) 

CALL GRID! 2,2) 

CALL CURVE! PHI ,Z, 359,0) 

CALL MESSAG! *THETA=25 DEGREE ' ,15, 11. ,8.5 ) 

CALL MESSAG! 'K=5. 5' ,5,11. ,8.2) 

CALL MESSAG! ' 01=1. 02=1. ' ,11,11. , 7. 9) 

CALL MESSAG! 'D3=l. D4=l. ' ,11 ,11 . ,7. 6 ) 

CALL MESSAG! ' P2=180 P3=0 P4=180 ' ,18,11. ,7.3 ) 
CALL MESSAG! ' A=l. B=l. C = l. ' ,14,11. ,7.0 ) 

CALL RESET( 'HEIGHT' ) 

CALL HEIGHT! .4) 

CALL MESSAG! 'PHI VS. DIRECTIVITY' ,19,3.1,9.8) 
CALL RESET! 'ALL' ) 

CALL NOCHEK 
CALL ENDGR(O) 

CALL PHYSOR! 11.85,7.6) 

CALL AREA20! 3 .5,2 .5 ) 

CALL FRAME 
CALL ENDPL(O) 

STOP 

END 



PROGRAM FOR THE CIRCULAR TWO-DIMENSIONAL GRAPH 



DIMENSION PHI! 360 ) ,R( 360 ), THETA! 360 ) ,Z( 360 ) ,X( 4 ) , Y( 4 ) 

REAL P4,P5>P6,P7,Wl,W2,W3,W4,W,G,H,K,Dl,D2,D3,D4,PH,TH,PI,a,b,c 
INTEGER P,P1,P2,P3,I,Q 



INPUT PARAMETERS 



PI=3. 1415917 

K=5.5 

Dl=l. 

D2=l. 

D3=l. 

D4=l. 

Pl= 180 
P2= 0 
P3= 180 
A = 1. 
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B = 1. 

C = 1. 

I = 25: USE P FOR P2DPCIR 



C 

C 

C 

C 

C 



c 



1 

200 

100 

C 

C — 
C 

c— 

c 



c 



MAIN PROGRAM 



P4=PI/180 

P5=P1*P4 

P6=P2*P4 

P7=P3*P4 

DO 100 1=0,360: MOVE C TO LINE NO. 174 FOR P2DICIR 
DO 200 P=0 >360 
TH=I*P4 
PH=P*P4 

W1=K*D1*SIN! PH )*COS! TH ) 

W2=K*D2*SIN! PH )*COS! TH )+P5 
W3=K*D3*SINt TH )*COS< PH )-P6 
W4=K*D4*SIN( TH )*COS< PH )+P7 
W=COS( W1 )+A*COS< W2 )+B*COS< H3 )+C*COS! H4 ) 
G=SIN( HI )-A*SIN( H2 )+B*SINl H3 )-C*SIN! W4 ) 
H=SQRT ( W*W+G*G )/( 1+A+B+C ) 

Z!P)=20*ALOG10(H) 

IF( Z( P ) . LE . -40 . ) THEN 
Z(P)=0. 

ELSE 

Z!P)=40.+Z(P) 

END IF 
THETA! P)=I 
PHI(P)=P 

HRITE ( 8,1 ) PHI(P),Z(P) 

FORMAT( 2X,2( E10 . 3 ,SX ) ) 

CONTINUE 

CONTINUE 



A PROGRAM FOR PLOTTING BY SHERPAOR BY TEK618 



CALL MEDBUF 
CALL COMPRS 

CALL TEK618: MOVE C TO THE UPPER LINE WHEN DRAWING WITH TEK618 
CALL PHYSOR( 1.5,1. ) 

CALL PAGE( 15.5,12. ) 

CALL AREA2D( 3. ,8. ) 

CALL COMPLX 
CALL HEIGHT( .2) 

CALL POLAR( 3.1415917/180. ,10. ,4. ,4. ) 

CALL GRID (1,1) 

CALL CURVE( PHI, Z, 360,0) 

CALL MESSAG( 'THETA=25 DEGREE ' ,15,9. 7,7. 7 ) 

CALL MESS AG ( 'K =5. 5' >5,9. 7,7. 4) 

CALL MESSAG( ' Dl=l . D2=l. 1 ,11 ,9. 7,7.1 ) 

CALL MESSAG( 'D3=l . D4=l . Ml , 9. 7 ,6 . 8 ) 

CALL MESSAG( ' P2=180 P3=0 P4=180 ' ,18,9. 7,6.5 ) 

CALL MESSAG!'A=1. B=l. C=1 . ' ,14, 9. 7,6 . 2 ) 

CALL MESSAG! 'O' ,1,3. 9,8.1 ) 

CALL MESSAG( '-10' ,3, 3. 7, 7.1) 

CALL MESSAG( '-20' ,3, 3. 7, 6.1) 

CALL MESSAG( ' -30 ' >3 >3 . 7 >5 . 1 ) 

CALL MESSAG! '-40 1 ,3,3. 7,4. 1 ) 

CALL MESSAG! ' 0( DEG )' ,6, 7. 6, 3. 9) 

CALL MESSAG( '180' ,3,0. ,3.9) 

CALL RESET( 'HEIGHT* ) 

CALL HEIGHT! .4) 

CALL MESSAG! 'PHI VS. DIRECTIVITY ' ,19,3 . 0 ,9 .5 ) 

CALL RESET! 'ALL' ) 

CALL NOCHEK 
CALL ENDGR(O) 

CALL PHYSOR! 10.95,6.8) 

CALL AREA2D! 3. 5, 2. 5 ) 

CALL FRAME 
CALL DONEPL 
STOP 
END 
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APPENDIX C 

PROGRAM FOR THE 3-D GRAPH USING SURMAT ABOUT X-Y PLANE 



xxxx x xxxxx xxxxx xxxxxxxx INFORMATION xxxx x xxxx x xxxxx x xxxx x xxxxxxx x xxxxx 

X X 

X THE OBJECTIVE OF THIS PROGRAM IS TO CALCULATE AND DRAW * 

x THE DIRECTIVITY PATTERN IN THREE-DIMENSION. x 

x THE DRAWING METHOD IS THE SUBROUTINE SURMAT OF DISSPLA. X 

x x 



X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

XXXXXXXX 



VARIABLE DEFINITION xxxxxxxxxxxxxxxxxxxxxxxxxxx 



X 

Z : DIRECTIVITY IN dB x 

ZMM : DOUBLE MATRIX TO GET X-Y PLANE x 

ZM : DUMMY VARIABLE OF Z x 

PI : PHASE (DEGREE) OF SOURCE A x 

P2 : PHASE (DEGREE) OF SOURCE B x 

P3 : PHASE (DEGREE) OF SOURCE C * 

P4 : VARIABLE TO CHANGE DEGREE TO RADIAN x 

P5 : RADIAL VALUE OF PI x 

P6 : RADIAL VALUE OF P2 X 

P7 : RADIAL VALUE OF PS x 

K : WAVE NUMBER x 

D1 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE 1 x 

D2 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE A x 

DS : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE B x 

D4 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE C x 

PHI : PHI IN DEGREE x 

I : THETA IN DEGREE x 

PH : PHI IN RADIAN x 

TH : THETA IN RADIAN x 

PI : 3 . 1415917 x 

x 

•xxxxxxxxxxxxx VARIABLE DECLARATION xxxxxxxx xxxxxxxxxxxxxxxxxx 
DIMENSION Z(81Q0),ZMM( 90,90) 

REAL P,T,R,X,Y,ZM,a,b,c 

REAL P4, P5 ,P6 ,P7, W1 , W2 ,W3 , W4,W,G,H,K,D1,D2,D3,D4, PH ,TH, PI 
INTEGER PHI,P1,P2,P3,I 



c 



c 



PROGRAM FOR THE UPPER HEMISPHERE 



c 

C 

c 

C 

C 

C 

C 



c 

c 

c 

c 

c 



A PROGRAM FOR PLOTTING BY SHERPAOR BY TEK618 



CALL COMPRS 

CALL TEK618: MOVE C TO THE UPPER LINE WHEN DRAWING WITH TEK618 
CALL RESET ( 3HALL ) 

CALL PAGE( 15. ,9. ) 

CALL PHYSORd. ,3. ) 

CALL SCMPLX 

CALL AREA2D(14.,8. ) 

CALL X3NAME ( 1H ,1) 

CALL Y3NAME ( 1H ,1) 

CALL Z3NAME ( 1H ,1) 

CALL VOLM3D (17., 17., 17. ) 

CALL VUANGL(+2Q. ,+20. ,500. ) 

CALL INTAXS 
CALL ZAXANG (-90.) 

CALL GRAF3D( 0. ,10. ,80. ,0. ,10. ,80. ,0. ,10. ,40. ) 

CALL BLSUR 

CALL BGNMAT( 90,90) 



INPUT PARAMETERS 



PI=3. 1415917 

K=PI 

Dl=l. 
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02 = 1 . 
DS=1. 
04=1. 
Pl= 180 
P2= 0 
P5= 180 
A = 1. 

8 = 0 . 

C = 0. 



C 

C 

C MAIN PROGRAM 

C 

C 

P4=PI/180 
P5=P 1*P4 
P6=P2*P4 
P7=P3*P4 
DO 200 I=-90,0 
TH=I*P4 

DO 100 PHI=1, 360,4 
PH=PHI*P4 

W1=K*D1*SIN! PH )*SIN( TH ) 

W2=K*D2*SIN( PH )*SIN( TH )-P5 
W3=K*D3*SINt TH )*COS( PH )-P6 
H4=K*D4*SINt TH )*COS( PH )+P7 
H=COS( HI )+A*COS( W2 )+B*COS( W3 )+C*COS( W4 ) 
G=SIN( HI ) -A*SINt W2 ) -B*SIN( W3 ) +C*SIN! W4 ) 
H=SQRT( W*W+G*G )/( 1+A+B+C ) 

Z( PHI )=20*ALOG10( H ) 

IF(Z( PHI KLE.-40. ) THEN 
Z( PHI )=0 . 

ELSE 

Z( PHI )=Z( PHI ) 

END IF 

R=40 . +Z( PHI ) 

P = FLOAT (PHI) 

T = FLOAT( I ) 

X=R*COS( TH )*COS( PH )+40. 

Y=R*COS! TH )*SIN! PH )+40 . 

ZM=-R*SIN( TH ) 

WRITE (8,1) X,Y,ZM 
1 FORMAT! 2X,3(E10.3,5X)) 

CALL GETMAT(X,Y,ZM,1,0) 

100 CONTINUE 

200 COfiTINUE 
C 

c 

C A PROGRAM FOR PLOTTING BY SHERPAOR BY TEK618 

C 

C 

CALL ENDMAT(ZMM,0) 

CALL SURMAT (ZMM, 1,90, 1,90,0) 

CALL DASH 
CALL MARKER! 17) 

CALL COMP LX 
CALL HEIGHT! .2) 

CALL SURVIS! ■BOTTOM 1 ) 

CALL MESSAG! 'DIRECTIVITY PATTERNS' ,19, 2.0, -.5) 
CALL RESET! 'COMPLX* ) 

CALL RESET! 'HEIGHT' ) 



c 

c 



PROGRAM FOR THE LOWER HEMISPHERE 



c 

C 

CALL BGNMAT! 90,90) 

C 

C INPUT PARAMETERS 

C - 

C 

PI=3. 1415917 

K=PI 

Dl=l. 

02 = 1 . 
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nnnnn 



MAIN PROGRAM 



D3=l. 
D<+=1. 
Pl= 180 
P2= 0 
P3= 180 
A = 1. 

B = 0. 

C = 0. 



P4=PI/180 
P5=P1*P4 
P6=P2*P4 
P7=P3*P4 

DO 250 I=-180,-90 
TH=I*P4 

DO 150 PHI=1, 360,4 
PH=PHI*P4 

W1=K*D1*SIN( PH )*C0S< TH ) 

W2=K*D2*SIN( PH )*COS! TH )+P5 
W3=K*D3*SIN! TH )*COS! PH )-P6 
W4=K*D4*SIN! TH )*COS( PH ) + P7 
W=COS( W1 )+A*COS( W2 )+B*C0S( H3 )+C*C0S! W4 ) 
G=SIN( W1 )-A*SIN! W2 )+B*SIN! H3 )-C*SIN! M ) 
H=SQRT( W*W+G*G )/( 1+A+B+C ) 

2! PHI )=20*AL0G10( H ) 

IF! 2! PHI ) . LE . -40. ) THEN 
2( PHI )=0 . 

ELSE 

2! PHI )=2! PHI ) 

END IF 

R=40 . +2! PHI ) 

P= FLOAT(PHI) 

T s FLOAT( I ) 

X=R*COS( TH )*COS( PH )+40 . 

Y=R*COS( TH )*SIN! PH )+40 . 

2M=R*SIN( TH ) 

WRITE! 8 >1 ) X,Y,2M 
1 FORMAT ( 2X>3( E10.3 ,5X ) ) 

CALL GETMAT!X,Y,2M,1>0) 

150 CONTINUE 

250 CONTINUE 

CALL ENDMAT! 2MM,0 ) 

CALL SURMAT ( 2MM, 1,90, 1,90,0 ) 

CALL DASH 
CALL MARKER! 17 ) 

CALL COMP LX 
CALL HEIGHT! .2) 

CALL SURVIS! ’BOTTOM' ) 

CALL ENDPLIO) 

CALL DONEPL 
CHECKSUM=012362 
STOP 
END 
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APPENDIX D 

PROGRAM FOR THE 3-DIMENSIONAL GRAPH USING VECTOR DRAWING MET 



x x * * * x* xxxx x x xx*x*xx x * * INFORMATION xx x xxx xxx xxx xx xxx x xxxxxxxxxxxx x x xx 

X X 

* THE OBJECTIVE OF THIS PROGRAM IS TO CALCULATE AND DRAW x 

* THE DIRECTIVITY PATTERN IN THREE -DIMENSION. * 

x THE DRAWING METHOD IS THE VECTOR DRAWING METHOD OF DISSPLA. * 

x x 



x x x xx x x xxxxxxx x xxxx x x x VARIABLE DEFINITION x x x xx xxxx x xx x xxxxxxxxxxx x x x 



X X 

x 2 : DIRECTIVITY IN dB * 
x XT : "TO" POINT OF X IN DRAWING VECTOR * 
x YT : "TO" POINT OF Y IN DRAWING VECTOR * 
x ZT : "TO” POINT OF Z IN DRAWING VECTOR x 
x XF : "FROM" POINT OF X IN DRAWING VECTOR x 
x YF : "FROM" POINT OF Y IN DRAWING VECTOR x 
x ZF : "FROM” POINT OF Z IN DRAWING VECTOR x 
x PI : PHASE (DEGREE) OF SOURCE A x 
x P2 : PHASE (DEGREE) OF SOURCE B x 
x P3 : PHASE (DEGREE) OF SOURCE C x 
x P4 : VARIABLE TO CHANGE DEGREE TO RADIAN x 
x P5 : RADIAL VALUE OF PI x 
x P6 : RADIAL VALUE. OF P2 x 
x P7 : RADIAL VALUE OF P3 x 
x K : WAVE NUMBER x 
x D1 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE 1 x 
x D2 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE A x 
x D3 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE B x 
x D4 : DISTANCE FROM THE GEOMETRICAL CENTER TO THE SOURCE C x 
x PHI : PHI IN DEGREE x 
x I : THETA IN DEGREE x 
x PH : PHI IN RADIAN x 
x TH : THETA IN RADIAN x 
x PI : 3. 1415917 x 
x x 



XXXXXXXXXXXXXXXXXXXXXX VARIABLE DECLARATION xxxxxxxxxxxxxxxxxxxxxxxxxx 
REAL R >Z >XT >YT >XF >YF >ZF >ZT >a >b >c 

REAL P4,P5»P6,P7»H1,H2,H3»H4>H,G,H»K>D1,D2»D3»D4»PH»TH»PX 
INTEGER PHI,P1,P2,P3,I 



DIMENSION Z( 0 : 2500 ) ,ZMM( 46 >46 ) , V( 2000 ) ,U( 2000 ) 

REAL P»T >R>X>Y >ZM >ZZM>A >B >C >V>U 

REAL P4,P5>P6,P7,H1,W2,W3,W4,W,G>H>K>D1>D2>D3>D4>PH>TH>PI 
INTEGER PHI>P1,P2,P3,I,S 



INPUT PARAMETERS 



PI=3. 1415917 
K=10 .55 
D1=0 . 5 
D2=0 .5 
D3=0 .5 
D4=0.5 
Pl= 0 
P2=180 
P3=180 
A = 1. 

B = 1. 

C = 1. 



MAIN PROGRAM 



P4=PI/180 
PS=P1XP4 
P6=P2xP4 
P7=P3xP4 
S = 1 
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DO 200 1=0,90 
TH=I*P4 

DO 100 PHI =0,180 
PH=PHI*P4 

W1=K*D1*SIN( PH )*SIN( TH ) 

- - W2=K*D2*SIN( PH )*SIN( TH )-P5 

- W3=K*D3*SIN( TH )*COS( PH )-P6 
W4=K*D4*SIN( TH )*COS( PH )+P7 

H=COS( HI )+A*COS( HZ )+B*COS( W3 )+C*COS( W4 ) 

G=SIN( HI )-A*SIN( W2 )-B*SIN( W3 HC*SIN( H4 ) 

H=SQRT( W*H+G*G )/( 1+A+B+C ) 

2( PHI )=20*ALOG10( H ) 

I F C ( Z( PHI ) . GT . -10 . ) . AND . ( Z( PHI ) . LT . -9 . ) ) THEN 
Z( PHI )= -30 . 

ELSE 

Z( PHI )= Z(PHI) 

END IF 
ZM=Z( PHI ) 

C WRITE ( 8,1 ) I ,PHI ,ZM 

IF (ZM .EQ. -30. ) THEN 
V(S) = FLOAT (PHI) 

U(S) = FLOAT (I) 

C WRITE (8,1) V( S ) ,U( S ) ,ZM 

S=S+1 
END IF 

C WRITE ( 8,1 ) I, PHI ,2 

Cl FORMAT( 2X,2( 13 ,SX ) ,F10 .5 ) 

C WRITE ( 8,1 ) AX( L ) ,BY ( L ) ,ZM 

Cl FORMAT ( 2X,3( F10.3,SX) ) 

100 CONTINUE 

200 CONTINUE 



A PROGRAM FOR PLOTTING BY SHERPA OR BY TEK618 



CALL MEDBUF 
C CALL COMPRS 
CALL TEK618 
CALL PHYSORd. ,1. ) 

CALL PAGE( 15.5,12. ) 

CALL AREA2D( 10.5,9.1) 

CALL FRAME 
CALL COMP LX 
CALL HEIGHT ( .2) 

CALL XNAME( ' PHI ( DEG ) $',10) 

CALL YNAME( *THETA( DEG )$ ' >10 ) 

CALL YTICKS( 2 ) 

CALL XTICKS( 9- ) 

CALL GRAF( 0,20 ,180 ,0,10,90) 

CALL CURVE ( V,U,S,-1 ) 

CALL MESSAG1 'THETA=25 DEGREE 1 ,15,11. ,8.5) 
CALL MESSAG( 'K=5.5 ' ,5,11 . >8. 2 ) 

CALL MESSAG( ' Dl=l. D2=l. 1 ,11,11. ,7. 9 ) 

CALL MESSAG( 'D3=l. D4=l. ',11,11. ,7. 6 ) 

CALL MESSAG( 'P2=180 P3=0 P4=180 * ,18,11. ,7.3 ) 
CALL MESSAG1 1 A=l. B=l. C=l. 1 ,14,11. ,7.0 ) 

CALL RESET ( 'HEIGHT* ) 

CALL HEIGHT ( .4) 

CALL MESSAG( 'TRACK( -20. DB) ' ,13,3.1,9.8) 

CALL RESET( ’ALL' ) 

CALL NOCHEK 
CALL ENDGR(O) 

CALL PHYSOR( 11.85,7.6) 

CALL AREA2D( 3 .5,2 .5 ) 

CALL FRAME 
CALL ENDPL(O) 

STOP 

END 
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nnnnn o onnno 



APPENDIX E 
P3D 



» x x x x»x«xxxxx*xx » x»xx m * INFORMATION xxx -* xxxx x* x ****x xxxxx x x»xxxx»xxxxx 

* 

OBJECTIVE OF THIS PROGRAM IS TO CALCULATE AND DRAW * 

DIRECTIVITY PATTERN IN THREE-DIMENSION. * 



* 

X 

X 

X 

* 

X 


THE 

THE 

THE 

ABOt 


* 




* 


2 : 


X 


ZMM 


* 


X : 


X 


Y : 


* 


ZM 


X 


PI 


X 


P2 


X 


P3 


* 


P4 


X 


P5 


X 


P6 


X 


P7 


* 


K : 


X 


D1 


* 


D2 


X 


D3 


* 


D4 


X 


PHI 


X 


I : 


* 


PH : 


* 


TH ; 


X 


PI : 


* 





DR 



m 



IG METHOD 
PLANE. 



IS THE SUBROUTINE SURMAT OF DISSPLA 



* 

* 

x 



VARIABLE DEFINITION 



DIRECTIVITY IN dB n . 

: DOUBLE MATRIX TO GET 0-<P plana 
DUMMY VARIABLE OF THETA 
DUMMY VARIABLE OF PHI 
DUMMY VARIABLE OF DIRECTIVITY 
PHASE (DEGREE) OF SOURCE A 
PHASE (DEGREE) OF SOURCE B 
PHASE (DEGREE) OF SOURCE C 
VARIABLE TO CHANGE DEGREE TO RADIAN 
RADIAL VALUE OF PI 
RADIAL VALUE OF P2 
RADIAL VALUE OF P3 
HAVE NUMBER 

DISTANCE FROM THE GEOMETRICAL 
DISTANCE FROM THE GEOMETRICAL 
DISTANCE FROM THE GEOMETRICAL 
DISTANCE FROM THE GEOMETRICAL 
: PHI IN DEGREE 
THETA IN DEGREE 
PHI IN RADIAN 
THETA IN RADIAN 
3.1415917 

VARIABLE DECLARATION 
DIMENSION Z( 2500 ),ZMM( 46,46) 

REAL X ,Y ,ZM >P4> P5 >P6 >P 7, HI ,H2 >H3 >H4,H,G,H,K,D1 ,D2 ,D3 >D4 , PH ,TH ,PI 
real a,b,c 

INTEGER PHI ,P1 >P2,P3 ,1 



CENTER TO THE SOURCE 1 
CENTER TO THE SOURCE A 
CENTER TO THE SOURCE B 
CENTER TO THE SOURCE C 



A PROGRAM FOR PLOTTING BY SHERPAOR BY TEK618 



CALL COMPRS 

CALL TEK618: MOVE C TO THE UPPER LINE HHEN DRAHING HITH TEK618 
CALL RESET (SHALL) 

CALL PHYSOR( 0.5,1. ) 

CALL PAGE( 15.5,12. ) 

CALL SCMPLX 

CALL AREA2D( 10.5,9.15) 

CALL X3NAME( 'THETAt DEGREE)' ,13) 

CALL Y3NAME( 'PHI( DEGREE)* ,11) 

CALL Z3NAME( 1 DIRECTIVITY ( DB ) ' > 15 ) 

CALL VOLM3D ( 17. ,17. ,17. ) 

CALL VUANGL( +20. ,+20. ,500. ) 

CALL INTAXS 
CALL ZAXANG (-90.) 

CALL GRAF3D( 0 . ,5. ,90. ,180. ,6. ,360. ,-40. ,2. ,0. ) 

CALL BLSUR 

CALL BGNMAT( 46,46) 



INPUT PARAMETERS 



PI=3. 1415917 
K=5 .5 
Dl=l. 

02 = 1 . 

D3=l. 

D4=l . 
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^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 



Pl = 180 
P2= 0 
P3= 180 
A * 1. 

B = 1. 

C = 1. 

C 

C 

C MAIN PROGRAM 

C 

C 

P4=PI/180 
P5=P1*P4 
P6=P2*P4 
P7=P3*P4 
DO 200 1=0,90,2 
TH=I*P4 

DO 100 PHI=180 ,360 ,4 
PH=PHI*P4 

N1=K*D1*SIN( PH )*SINt TH ) 

W2=K*D2*SIN( PH )*SIN(TH)-P5 

N3=K*D3*SIN( TH )*COSt PH )-P6 

H4=K*D4*SIN( TH )*COSt PH )+P7 

H=COS( HI )+A*COS( N2 )+B*COS( N3 )+C*COS( W4 ) 

G=SIN( HI )-A*SIN( W2 )-B*SIN( H3 ) +C*SINt W4 ) 

H=SQRT ( W*W+G#G )/( 1+A+B+C ) 

Z(PHI)=20*ALOG10(H) 

I F C Z ( PHI ) . LE . -40 . ) THEN 
Z(PHI)=-40. 

ELSE 

Z( PHI )= Z(PHI) 

END IF 
X=I 
Y=PHI 
ZM=Z( PHI ) 

WRITE (8,1) X,Y,ZM 
1 FORMATt 2X, 3 (E10. 3, 5X)) 

CALL GETMAT(X,Y,ZM,1,0) 

100 CONTINUE 

200 CONTINUE 
C 

C 

C A PROGRAM FOR PLOTTING BY SHERPAOR BY TEK618 

C 

C 

CALL ENDMATt ZMM ,0 ) 

CALL SURMAT ( ZMM,1 ,46 ,1,46 ,0 ) 

CALL DASH 
CALL MARKER( 17 ) 

CALL COMPLX 
CALL HEIGHT ( .2) 

CALL SURVISt •BOTTOM 1 ) 

CALL MESSAG( ’K=5.5’ ,5,11. ,8.5) 

CALL MESSAGt 'Dl=l. D2=l. ' ,11,11. ,8. 2 ) 

CALL MESSAGt 'D3=l. D4=l . ' ,11 ,11 . , 7 . 9 ) 

CALL MESSAGt 'VIEW ANGLE= 20, 20 ' ,18 ,11 . ,7.6 ) 

CALL MESSAGt 'P2=180 P3=0 P4=180 ' ,18 ,11 . ,7.3 ) 

CALL MESSAGt ‘A=l. B=l. C=l. ' ,14,11. ,7.0 ) 

CALL RESETt ‘HEIGHT’ ) 

CALL HEIGHTt .4 ) 

CALL MESSAGt 'DIRECTIVITY PATTERN',19, 3. 05,9.8 ) 
CALL ENDGRtO) 

CALL PHYSORt 11.35,7.6) 

CALL AREA2D( 3 .5 ,2 .5 ) 

CALL FRAME 
CALL DONEPL 
CHECKSUM=012362 
STOP 
END 
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APPENDIX F 
TRACE 



xxxxxxxx x x xxxx xxkxxxxxx INFORMATION xxxxxxxxxx x xxxxxx x xxxxxxx: 

X 

x THE OBJECTIVE OF THIS PROGRAM IS TO CALCULATE AND DRAW 

x THE CONSTANT dB LEVEL OF THE THREE-DIMENSIONAL PLOT, 

x THE DRAWING METHOD IS THE SUBROUTINE CURVE OF DISSPLA. 

x 

xxxxx x xxxxx xxxxxxxxxxx VARIABLE DEFINITION x x x x xxxxx xxx xx xx x x xxxx x xx xx 



X 
















X 


X 


Z : 


DIRECTIVITY IN dB 












X 


X 


v : 


ARRAY TO TRACK 3> 












X 


X 


u : 


ARRAY TO TRACK U 












X 


X 


PI 


PHASE (DEGREE) OF SOURCE A 












X 


X 


P2 


PHASE (DEGREE) OF SOURCE B 












X 


X 


P3 


PHASE (DEGREE) OF SOURCE C 












X 


X 


P4 


VARIABLE TO CHANGE DEGREE TO 1 


RADIAN 










X 


X 


P5 


RADIAL VALUE OF PI 












X 


X 


P6 


RADIAL VALUE OF P2 












X 


X 


P7 


RADIAL VALUE OF P3 












X 


X 


K : 


WAVE NUMBER 












X 


X 


01 : 


: DISTANCE FROM THE GEOMETRICAL 


CENTER 


TO 


THE 


SOURCE 


1 


X 


X 


D2 : 


: DISTANCE FROM THE GEOMETRICAL 


CENTER 


TO 


THE 


SOURCE 


A 


X 


X 


D3 j 


: DISTANCE' FROM THE GEOMETRICAL 


CENTER 


TO 


THE 


SOURCE 


B 


X 


X 


D4 ; 


; DISTANCE FROM THE GEOMETRICAL 


CENTER 


TO 


THE 


SOURCE 


C 


X 


X 


PHI 


: PHI IN DEGREE 












X 


X 


I : 


THETA IN DEGREE 












X 


X 


PH : 


: PHI IN RADIAN 












X 


X 


TH : 


; THETA IN RADIAN 












X 


X 


PI : 


: 3.1415917 












X 


X 


xxxxxxx* 


VARIABLE DECLARATIOI 


M xxxx* 










X 



DIMENSION Z( 0 : 2500 ) >ZMM( 46 ,46 ) , V( 2000 ) ,U( 2000 ) 

REAL P,T,R,X,Y,ZM,ZZM,A,B,C,V,U 

REAL P4,P5,P6,P7,W1,W2,W3,W4,W,G,H,K,D1,D2,D3,D4,PH,TH,PI 
INTEGER PHI,P1,P2>P3>I>S 
c 

C 

C INPUT PARAMETERS 

C 

C 

PI=3. 1415917 
K=10 .55 
Dl=0.5 
02=0.5 
03=0.5 
D4=0 .5 
Pl= 0 
P2=180 
P3=180 
A = 1. 

B = 1. 

C = 1. 

C 

C 

C MAIN PROGRAM 

C 

C 

P4=PI/180 

P5=P1*P4 

P6=P2XP4 

P7=P3*P4 

S=1 

DO 200 1=0,90 
TH=IXP4 

DO 100 PHI=0 , 180 
PH=PHI*P4 

W1=K*D1*SIN( PH )*SIN( TH ) 

W2=KXD2XSIN( PH )xSIN( TH )-P5 
W3=K*D3*SIN( TH )xCOS( PH )-P6 
W4=KXD4XSIN( TH )xCOS( PH )+P7 
W=COS( W1 )+A*COS( W2 )+BxCOS( W3 )+C*COS( W4 ) 
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* * * * * 



ooooo 



G=SIN( W1 )-A*SIN( W2 )-B*SIN< W3 )+C*SIN( W4 ) 

H=SQRT( W*W+G*G )/( 1+A+B+C ) 

Z(PHI)=20*AL0G10(H) 

IF( ( Z( PHI ) . GT . -10 . ) . AND . ( Z( PHI ) . LT . -9 . ) ) THEN 
Z( PHI )= -30. 

ELSE 

Z( PHI ) s Z(PHI) 

END IF 
ZM=Z( PHI ) 

C WRITE (8,1) I >PHI >ZM 

IF (ZM .EQ. -30. ) THEN 
V(S) = FLOAT (PHI) 

U(S) = FLOAT (I) 

C WRITE ( 8,1 ) V(S),U(S),ZM 

S=S+1 
END IF 

C WRITE ( 8,1 ) I,PHI,Z 

Cl FORMAT ( 2X,2( I3,5X) ,F10.5) 

C WRITE ( 8,1 ) AX( L ) >BY( L ) ,ZM 

Cl FORMAT ( 2X,3( F10.3,5X) ) 

100 CONTINUE 

200 CONTINUE 



A PROGRAM FOR PLOTTING BY SHERPA OR BY TEK618 



CALL MEDBUF 

C CALL COMPRS>MOVE C TO TEK619 TO PRINT USING SHERPA 
CALL TEK618 
CALL PHYSOR( 1. >1. ) 

CALL PAGE( 15.5,12. ) 

CALL AREA2Df 10.5^9.1) 

CALL FRAME 
CALL COMP LX 
CALL HEIGHT ( .2) 

CALL XNAME( 'PHI(DEG) $',10) 

CALL YNAME( 'THETA( DEG )$ ' ,10 ) 

CALL YTICKS( 2 ) 

CALL XTICKS(^) 

CALL GRAF( 0,20,180,0,10, 90 ) 

CALL CURVE (V,U,S,-1) 

CALL MESSAG( 'THETA=25 DEGREE * ,15,11. ,8.5) 

CALL MESSAG( 'K=5.5 1 ,5,11 . ,8. 2 ) 

CALL MESSAG( 'Dl=l. D2=l. ‘ ,11,11. ,7. 9) 

CALL MESSAG( *D3=1. D^f=l. 1 ,11,11. ,7.6 ) 

CALL MESSAG( 'P2=180 P3=0 P^=180 ' ,18,11. ,7.3 ) 

CALL MESSAG( 'A=l. B=l. C=l. 1 ,14,11. ,7.0 ) 

CALL RESET ( 'HEIGHT' ) 

CALL HEIGHT ( .4) 

CALL MESSAGt 'TRACK* -20. DB) ' ,13,3.1,9.8) 

CALL RESET ( 'ALL' ) 

CALL NOCHEK 
CALL ENDGR(O) 

CALL PHYSOR( 11.85,7.6) 

CALL AREA2D( 3. 5, 2. 5) 

CALL FRAME 
CALL ENDPL(O) 

STOP 

END 
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APPENDIX G 

LIST OF COMPUTER OUTPUT 



• SAMPLE CASE 1 



EXXEEXXEE X EE E E X XE EX XEEXXXXXEXXXXXEXEX XEEXEEE X XEX XXXXEXXEXXXXXXXXXXXXKXXXXXXEXXXXXXXXXXXXXXXXXXX 

A B C K 0, D-, D, D. (D, <0, <p, <P (3 H F 

X EEX X X EXEXXXXXXXEXXXEXXXXXXXXXXXEEXXKXXEEXTXXXXEEXE X XKEXXXXXXXEXXXXXXXXXEEXXXEXXXXEEXXXXEXXXXXE 



1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


31 


0.955567 


- 0.394776 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


32 


0.970322 


- 0.261679 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


33 


0.982069 


- 0.157164 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


34 


0.990826 


- 0.080049 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


35 


0.996641 


- 0.029226 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


36 


0.999579 


- 0.003656 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


37 


0.999728 


- 0.002363 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


38 


0.997193 


- 0.024419 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


39 


0.992094 


- 0.068947 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


40 


0.984566 


- 0.135104 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


41 


0.974756 


- 0.222079 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


42 


0.962821 


- 0.329093 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


41 


0.948922 


- 0.455388 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


44 


0.933230 


- 0.600225 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


45 


0.915917 


- 0.762881 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


46 


0.897156 


- 0.942645 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


47 


0.877121 


- 1.138807 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


48 


0.855985 


- 1.350672 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


49 


0.833918 


- 1.577533 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


50 


0.811083 


- 1.818697 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


81 


0.261312 


- 11.656804 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


82 


0 . 255355 


- 11.857100 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


83 


0.250133 


- 12.036570 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


84 


0.245633 


- 12.194245 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


85 


0 . 241845 


- 12.329250 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


86 


0.238759 


- 12.403822 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


87 


0.236367 


- 12.528270 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


88 


0.234663 


- 12.591125 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


89 


0.233642 


- 12.628973 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


90 


0.233302 


- 12.641622 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


91 


0.233642 


- 12.628973 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


92 


0.234663 


- 12.591125 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


93 


0.236367 


- 12.528270 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


94 


0.238759 


- 12.440822 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


95 


0.241845 


- 12.329267 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


96 


0.245633 


- 12.194263 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


97 


0.250133 


- 12.036578 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


98 


0 . 255355 


- 11.857116 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


99 


0.261311 


- 11.656837 


1.0 


1.0 


1.0 


10.55 


0.50 


0.50 


0.50 


0.50 


0 


180 


180 


0 


100 


0.268014 


- 11.436835 
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